Generalized Inverse Participation Ratio as a Possible Measure of Localization for 

Interacting Systems 
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We test the usefulness of a generalized inverse participation ratio (GIPR) as a measure of Anderson 
localization. The GIPR differs from the usual inverse participation ratio in that it is constructed 
from the local density of states rather than the single-electron wavefunctions. This makes it suitable 
for application to many-body systems. We benchmark the GIPR by performing a finite-size scaling 
analysis of a disordered, noninteracting, three-dimensional tight-binding lattice. We find values for 
the critical disorder and critical exponents that are in agreement with published values. 
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I. INTRODUCTION 

Anderson localization is a phenomenon in which quan- 
tum particles may be localized due to a random potential, 
even though the particles are classically unbound. 1 The 
theory for noninteracting particles is well-developed: in 
one and two dimensions, particles are localized by ar- 
bitrarily weak disorder, and in three dimensions states 
may be localized or extended depending on the strength 
of disorder^ 

Most real particles are interacting, however, and there 
has been an ongoing effort to understand how interac- 
tions modify the noninteracting picture, either because 
of screening of the disorder potential, or because of loss 
of quantum coherence due to inelastic scattering. 3 Until 
recently, neither of these effects was believed sufficient to 
change the noninteracting picture at zero temperature. 
However, experiments 4 - in two dimensional semiconduc- 
tor films identified a zero-temperature metal-insulator 
transition (MIT) that appears to result from electron 
interactions £ More recently, it has been suggested that 
weakly-interacting one and two-dimensional systems will 
exhibit a finite-T Anderson MIT&L 

There is also interest in Anderson localization in 
strongly interacting systems^ Many of the most interest- 
ing strongly correlated materials are insulators, but can 
have their electronic properties tuned by chemical dop- 
ing. Of particular interest are materials, such as the high 
temperature superconductors, whose parent compounds 
have an interaction-driven Mott insulating phase. These 
materials become superconductors when doped with a 
few percent of electron or hole donor atoms, but pass 
through various intermediate phases in which disorder 
seems to play an important role. There is an abundancy 
of questions about how the electronic properties of these 
materials are modified by doping-related disorder. Of 
particular relevance to this work, there have been recent 
questions about how localization physics is altered near 
the Mott MIT,Sr— and about the phase transition be- 
tween the Anderson and Mott insulating phases* 1 ^— 

Finally, trapped atomic gases in random optical lat- 



tices have now been experimentally realized^ - — These 
systems are interesting because the strength of the atom- 
atom interactions can by tuned by application of an ex- 
ternal magnetic field. There is therefore the prospect of 
making a controlled study of Anderson localization as a 
function of interaction strength. 

Numerical calculations have played an important role 
in understanding Anderson localization in noninteract- 
ing systems. However, many of the techniques devel- 
oped for measuring localization in noninteracting sys- 
tems cannot be extended to interacting systems since 
they require knowledge of the single-particle eigenstates 
of the system and, with the exception of self-consistent 
field calculations, many-body wavefunctions cannot gen- 
erally be written as a simple product of single-particle 
states. There is, therefore, an interest in developing new 
numerical methods for studying the Anderson MIT in 
interacting systems. 

With this in mind, there have been several proposals 
that the localization transition can be detected by study- 
ing the statistical properties of the local density of states 
(LDOS) p(r, lj). The geometric average of the LDOS, 
p g (co), is an order parameter for the Anderson MIT in the 
limit of infinite system siz o 23 ' 24 because it vanishes when 
the local spectrum is discrete. In infinite systems, this oc- 
curs only at energies at which the states are localized and 
not at which the states are extended. A generalization 
of dynamical mean field theory based on incorporating 
p g {oj) into the self-consistency cycle was developed to 
study interacting disordered systems j 15 ' 16 ' 24 As a prac- 
tical measure of localization in finite systems, however, 
p g {oj) is problematic because the spectrum is always dis- 
crete, and this can obscure the Anderson MIT^^ More 
recently, several groups have suggested that the Anderson 
transition can be detected by studying the distribution of 
p(r, uj) values , 16 ' 27 " — and it has been shown that this dis- 
tribution scales differently with system size for localized 
and delocalized states^ 

In this work, we consider a quantity, the generalized 
inverse participation ratio (GIPR), that is related to the 
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LDOS via 
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Equation (P) is denned for a lattice, so that p(r i; w) is 
the density of states projected onto the local Wannier 
orbital at the ith site of the lattice. The GIPR was used 
previously in finite size scaling studies^ but a careful 
examination of its scaling properties has not been made. 
This is the purpose of this paper. 

The GIPR is analogous to the usual inverse participa- 
tion ratio (IPR) for noninteracting systems, 
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where ^ a (ri) is a single-particle wavefunction with quan- 
tum number a in the basis of Wannier orbitals. The IPR 
is conventionally defined with q — 2 and can be used to 
distinguish Anderson localized and extended states: for a 
finite d-dimensional system of linear size L, I<i, a satisfies 



lim h. a 



_ J 1/L d (extended states) 
const, (localized states), 



(3) 



for states that are far from the Anderson MIT, and ex- 
hibits multifractal scaling ) 23 ! 30 " — 



lim 7 2 , Q = L- d *F[(W - WJL 1 '"] 



(4) 



near the transition. Here, d,2 is the fractal dimension for 
q = 2, v is a critical exponent, and W and W c are the 
disorder and critical disorder strengths respectively. 

For noninteracting systems, G2(u>) reduces to the IPR 
when u> is equal to one of the eigenenergies of the system. 
This follows from substituting 



p(r i ,cj)=J2\* a (r i )\ 2 S(u;-E a ), 



(5) 



into Eq. (fTJ), where E a are the discrete eigenenergies of 
the disordered lattice. However, for a general value of uj 
not equal to one of the eigenenergies, G-z(oS) is not well 
defined if the 5-functions in Eq. ((5]) are infinitely sharp, 
and the relationship between the IPR and the GIPR is 
therefore ambiguous. Moreover, we show below that if 
one broadens the ^-functions by an amount 7, there is no 
limiting value of 7 in which the GIPR reduces to the IPR. 
The goal of this paper is to demonstrate that the GIPR 
can nonetheless be used to detect the Anderson MIT and 
to determine the critical parameters W c , d<x and v. 

We benchmark the GIPR by performing finite size scal- 
ing for a disordered noninteracting model, where the crit- 
ical properties are well known. In Sec. [HI we discuss how 
the broadening of the <5-functions in Eq. ([5]) is expected 
to affect the finite size scaling, and use this to select an 
optimal broadening. In Sec. IIII1 we show the results of 
numerical finite size scaling, from which we extract values 
for the critical disorder and critical exponents at the An- 
derson MIT. We show that, with an appropriate choice 
for 7, it is possible to extract critical properties. 



II. CALCULATIONS 



The noninteracting Anderson model is 



(6) 
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where \i) is the ket for a Wannier orbital at position i 
on the lattice and indicates that the sum is over 

nearest-neighbour sites. The hopping matrix element is 
taken to be t = 1, and it therefore sets the energy scale, 
while the site energies are taken from a uniform dis- 
tribution of random values ranging from — W/2 to W/2, 
where W is the strength of disorder. Calculations are 
performed for a three-dimensional (d — 3) cubic lattice 
of linear size L and with N s = L 3 lattice points. 

We use a recursion method 3 - 3 - to find the local Green's 
function G(r,-, oj+ij) at site i, where 7 is a small but finite 
shift off the real frequency axis. This method introduces 
an error through truncation of the recursion algorithm, 
and we have been careful to adjust the truncation crite- 
rion so that this error is much smaller than the error due 
to disorder averaging. The LDOS is given by the imag- 
inary part of G(ri,ui + £7). Formally, this is equivalent 
to 



p 7 (r 
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where E a are the eigenenergies for a particular disorder 
realisation. Once /? 7 (rj,o;) is known, the GIPR is calcu- 
lated from Eq. (TT]). In this work, we focus on the band 
center (uj = 0), where the Anderson transition is well- 
characterized. In particular, the Anderson MIT occurs 
at a critical disorder W c — 16. 5t for the uniform disorder 
distribution used herei 34 ' 35 

One of the main issues we face is how to choose 7. In 
the remainder of this section, we discuss how this choice 
affects both the LDOS and the scaling behavior of the 
GIPR. The relevant energy scale for comparison is the 
level spacing at the band center, A = 1/pqN s , where 
po is the system-averaged density of states at w = 0. 
For strongly disordered systems, A w W/N s , while for 
weakly disordered systems, A w D/N s , where D is the 
bandwidth of the disorder-free lattice. For the cubic lat- 
tice considered here, the Anderson transition occurs at 
an intermediate disorder strength, so that A lies between 
these two limits. 

Figure [1] shows the dependence of the LDOS on 7. 
When 7 > A, the LDOS at oj is an average over states 
with \oj — E a \ < 7, with the consequence that the LDOS 
is more spatially uniform than the individual eigenstates 
making up the LDOS. Thus in Fig. [IJ the sites A and 
B are spatially separated, and both have broad peaks at 
uj = 1. It is not possible to tell, based on the LDOS 
for 7 = W/N S) whether these peaks indicate a single 
eigenstate or a cluster of eigenstates that happen to be 
close in energy. It is only when 7 A that we see 
that the local spectrum is quite different at the two sites 
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FIG. 1: (Color online) Local density of states for two well- 
separated lattice sites, "A" and "B", in a disordered lattice. 
Panels show the LDOS at (a)-(c) A and (d)-(f) B. All spectra 
are for the same configuration of disorder, but have different 
values of 7. Results are for (a), (d) 7 = 0.01W/N S ; (b), (e) 
7 = W/N s ; (c), (f) 7 = 2W/N S , where the lattice has iV s = 4 3 
sites and W = 13. Insets show expanded views of the LDOS 
near uj = 1. 



[Figs. HJa) and (d)]. This suggests that the finite size 
scaling of the LDOS, and by extension the GIPR, should 
do a better job of distinguishing localized and extended 
states as 7 is reduced. 

However, the fact that the spectrum develops discrete 
peaks when 7 < A does not mean that the LDOS sam- 
ples only individual eigenstates in this limit. This is be- 
cause most energies do not coincide with a peak. When 
7 <§; A, the energy uj = lies in the tails of the surround- 
ing peaks and Eq. J7| becomes 



p 7 (r l7 0) = ^^ 

TT * ■ 



E 2 



(8) 



This means that even in the limit 7—^0, /?(rj,0) is av- 
eraged over a nonzero number of states. The LDOS at 
u) = 0, and by extension the GIPR, does not change 
qualitatively when 7 is reduced much below A. 

We can learn more about the GIPR scaling by substi- 
tuting Eq. ([7]) into Eq. ([1} , from which we obtain 



where 



{El 



(9) 



(10) 



is a weighting factor satsifying w a — 1. In the limit 
of vanishing disorder, the wave functions are plane waves 
with l^a^)! 2 = N-\ and Eq. © gives G 2 (0) = AT" 1 ; 



this result is independent of 7 and is identical to the 
scaling result for the IPR. 

In the limit of large disorder, W 3> W c , it is useful to 
rearrange Eq. @ to obtain, 

G 2 (0) = Y, w lh a + Y, w ^Y,^^)\ 2 \^f3(r t )\ 2 . 

a 0,^,9 i 

(11) 

The first term on the right hand side is a weighted sum of 
IPR values for eigenstates with \E a \ < 7, while the sec- 
ond term consists of cross terms between pairs of eigen- 
states. The second term can be neglected when the dis- 
tances between these localized states are large compared 
to the localization length £. We can estimate the typical 
distance between centers of localization of the states in 
Eq. (jTTJ) for the case 7 > A. In this case, there are of or- 
der 27/ A states with \E a \ < 7, and the mean separation 
of these is 



£~ L(A/2 7 ) 1 / d . 



(12) 



The product |^ Q (rj)| 2 |^' ( g(rj)| 2 for two states separated 
by £ has a maximal value of order exp[— 2L(A/2 r y) 1 ' d /£], 
at the midpoint between the centers of localization. It 
follows that the second term in Eq. (fTTj) vanishes for 
L/£ — > 00, in which limit the GIPR is expected to scale 
like the IPR. 

For finite L, however, the second term in Eq. (1111) in- 
troduces finite size corrections to the GIPR that make it 
scale differently from the conventional IPR. In order to 
minimize these corrections, we want to make £ as large as 
possible, which is achieved by taking 7 as small as possi- 
ble. We emphasize, however, that Eq. ([12")) only holds for 
7 > A, and that w a is independent of 7 when 7 -C A, 
namely 



lim w a 

7^0 



E- 



(13) 



In other words, I ceases to increase when 7 is much less 
than A. Our analysis therefore suggests that one cannot 
do much better at minimizing finite size effects than by 
taking 7 ~ A. 

Finally, having established that 62(0) is determined 
by the first term in Eq. (fTT|) when L/£ 3> 1, we show 
that the weighting terms do not affect the GIPR scaling 
in this limit. We write I% ia « It{E a ), where l2{E) is a 
slowly varying function of E near E = 0, so that 



G a (0)«J a (0)XN 



(14) 



For 7 > A, we may estimate the sum over eigenstates by 
A- 1 J dE(E 2 + j 2 )- 2 



E«4 



[A- 1 $dE(E 2 +7 2 )- 1 ] 
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(15) 



Taking 7 oc iV s 1 eliminates the L-dependence of the 
weighting factors in Eq. (TTT1) . 
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FIG. 2: (Color online) Probability distribution of the loga- 
rithm of the local density of states at uj — 0. The effect of 7 
on Pp(hi p 7 ) is shown for (a) extended and (b) localized states 
for a fixed system size. The effect of system size on P P (ln p 7 ) is 
shown for (c) extended and (d) localized states for 7 = W/N s . 
Here, p 7 is the normalized LDOS, p 7 = p 7 /(p 7 ), where (p 7 ) 
is the system-averaged LDOS at u — 0. Results are shown 
for 16 (L = W),5(L= 15), and 2 (L = 20) disorder configu- 
rations, such that the number of LDOS values in each case is 
roughly the same. 



In summary, we have shown that the GIPR will repro- 
duce the scaling of the IPR in the limits of vanishingly 
weak and strong disorder. Moreover, we have shown that 
using smaller 7 values to calculate the GIPR is prefer- 
able, down to 7 ~ A. Many numerical methods converge 
faster for larger 7 and, for these, 7 ~ A will be optimal. 
In the next section, we examine whether the finite size 
effects near W c limit our ability to extract the critical 
behaviour. 




FIG. 3: (Color online) Effect of 7 on the probability distri- 
bution of the GIPR at uj = 0. Results are for (a) extended 
(W = 13) and (b) localized (W = 20) states, and are for 1500 
disorder configurations with L — 10. 



almost independent of L for the extended case (W = 
13). The small leftward shift seen in the extended case 
is presumably due to finite-size effects, which are more 
pronounced here than in Ref. HU. Despite its smallness, 
this leftward shift is problematic because it obscures the 
signature of the Anderson transition in P p (lnp 7 ). This is 
a potentially important issue for many-body calculations 
where accessible system sizes tend to be severely limited. 
It appears that, as with other measures of localization, 
the usefulness of the LDOS distribution will depend on 
the inclusion of finite-size corrections. 

Figure [3] shows the probability distribution function 
Pc^lnG^) for the GIPR, obtained by calculating G2(u>) 
at lu = for 1500 distinct impurity configurations. This 
figure shows that the width of the distribution depends 
strongly on 7, and that Pb(l n G2) is sharply peaked when 
7 > A. Because the distribution of lnG2 is narrow, the 
mean and most probable values of the distribution are 
close to each other. For this reason we study the finite 
size scaling of the typical GIPRj2£ 



G| yp (a;) = exp[<mG 2 (a;))] 



(16) 



III. RESULTS 

We plot, in Fig. [5J the probability distribution of the 
logarithm of the LDOS at ui = for different values of 
7 and for different system sizes. Figures [Ha) and (b) 
show that 7 affects both the peak position and shape of 
the distribution. In particular, the peak position of the 
distribution P p (lnp 7 ) is proportional to 7 for 7 < W/N s , 
in accordance with Eq. ((5J). For 7 > W/N s , the peak 
position and width are weak functions of 7. 

In Figs.^c) and (d), we show the L-dependence of the 
distribution of the normalized LDOS, p 1 = p 7 /(p 7 ) with 
(p-y) the sample-averaged density of states. Schubert et 
ali2£ showed that the scaling of the distribution P p (ln p 7 ) 
can be used to distinguish localized and extended states: 
P p (In p 7 ) shifts to the left with increasing L for localized 
states, and is independent of L for extended states. Here, 
we find that there is indeed a pronounced shift for the 
localized case (W = 20), and that the distribution is 



where (...) refers to an average over disorder configura- 
tions. 

We argued in the previous section that one should take 
7 oc A, where A depends on both N s and W. To under- 
stand whether the FT-dependence of A is important, we 
take two cases: 7 oc W/N s and 7 oc W c /N s , where W c 
here refers to the accepted value of 16.5. As we dis- 
cussed in Sec. HH the first choice overestimates the In- 
dependence of A, while the second choice underestimates 
it. Note that there is nothing fundamental about the pro- 
portionality constant W c in the second case; it was chosen 
because it gives 7 values that are quantitatively close to 
those in the first case. In total, we have taken four cases: 
two with 7 oc W/N s (7 = W/N s and 7 = 2W/N S ) and 
two with 7 oc W c /N s (7 = W c /N s , and 7 = 2W C /N S ). 

Figure [4ja) shows the dependence of G2 yp (0) on L for 
different strengths of disorder for the case 7 = W/N s . 
At short length scales, all the systems are in the criti- 
cal region (albeit in a region where finite size corrections 
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FIG. 4: (Color online) Scaling of the GIPR for 7 = W/N s . 
(a) Plots of G 2 yP as a function of L for 3000 disorder config- 



urations, (b) Y2, denned by Eq. (|18[) . for best fit values of 
Ao, d,2 and y. For these parameters, the critical disorder at 
which all curves cross is W c = 16.2, corresponding to the op- 
timized fitting parameters Ao — 0.64, y = 1.6, and d2 = 1.3. 
(c) Plot of dYn/dW at W — W c (symbols), along with power 
law fits (solid lines) to the data. The exponents in the fit- 
ted curves give v~ 1 , from Eq. (|2U|) . The two outlying curves 
are for the extremal values (d,2,W c ,y) = (1.5,15.9,3.2) and 
(1.2, 16.4, 1.2) and are used to determine uncertainties for v. 
The middle curve is for the optimized parameters, from which 
we obtain v = 1.1. Results are summarized in Table U 



are significant), and therefore all show similar size depen- 
dence. At long length scales, however, the lines diverge. 
For W < W c , the slope becomes steeper with increasing 
L, consistent with a crossover to L d with d = 3. For 
W > W c , the slope decreases with increasing L, consis- 
tent with a crossover to a constant value. This figure 
suggests that the GIPR is indeed able to distinguish lo- 
calized and extended states, even for the relatively small 
systems studied here. 

We show that the GIPR displays the same critical be- 
haviour as the IPR near W c , namely that 
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typ _ T -d,2 
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(W - W c )L l l v 



Ao 
Ly 



L 2 y 



(17) 

where Aj are finite size corrections and y is the critical 
exponent for the leading-order irrelevant variable^ In 
all cases, we are able to obtain good scaling behaviour 
for 4 < L < 17 with Aj = for j > 1. We thus have five 
fitting parameters: d 2 , W c , u, A , and y. 

We now describe the fitting procedure, using the case 
7 = W/N s as an example. Figure B|b) shows a plot of 



G 2 yp L d 



A, 

Ly 



(18) 



versus W for the optimal values of Ao, y, and d 2 . Error 
bars on the data are the root-mean-square uncertainty 
in due to the finite width of the GIPR distributions 
(shown, e.g., in Fig. [3]). The solid curves in Fig. IDJb) 
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FIG. 5: (Color online) Critical parameters for (a) 7 = W/N s 
and (b) 7 = W c /N s at the best-fit values (a) y — 1.6 and (b) 
y — 2.9. Intensity scale shows the value of Wx that minimizes 
X 2 locally for each and Ao. Circles indicate best-fit values 
of d2 and Ao, obtained from the global minimum ofx 2 - Black 
contours bound the regions Xrcd < !• 



are cubic fits to the data points. Each pair of curves 
crosses at a different disorder strength, denoted Wj ± 
SWj, where j e [1, A^ cross ] and N cmss is the the number of 
such crossing points. (For the 8 curves shown in Fig.QJb), 
there are N cross = 28 crossing points.) The uncertainties 
SWj are calculated from the uncertainties in the fitting 
parameters. If the scaling form Eq. (fl7j) holds and the 
critical parameters are correctly chosen, all curves will 
cross at a single point, Wx- For each Ao, y and d 2 , we 
find Wx(Ao,d,2,y) by minimizing 



JV crraB 

E 



Wj - W x 
SWi 



(19) 



Plots of Wx(A ,d 27 y) are shown in Fig. [5] for optimal 
values of y for 7 = W/N s and for 7 = W c /N s . We 
extract our own best-fit values for W c from the global 
minima of % 2 (Ao, d 2 , y), and these are shown as circles 
in Fig. [SJ A qualitative sense of the goodness-of-fit can 
be obtained from Fig. 0|b), which is based on the best- 
fit parameters for 7 = W/N s . A quantitative measure 
of goodness-of-fit can be obtained from the reduced chi- 
square vjr cd = x 2 /(A r c — 1). Figure [5] shows contours 
around the region of parameter space Xrcd < 1- I n this 
region, all ^(^0 curves cross, within error, at a common 
point. 

The best-fit values for W c and d 2 are summarized in 
Table HI along with previously published values. Quanti- 
ties in brackets are extremal parameter values satisfying 
Xrcd < 1; an d are used to estimate the uncertainty in 
the critical parameters. The values for W c and d 2 found 
from this analysis are generally within uncertainty of the 
previously published results. 
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7 


d 2 


W c 


V 


y 


W/N 3 


1.3 (1.2,1.5) 


16.2 (15.9,16.4) 


1.1 (1.0,1.1) 


(1.2,3.2) 


2W/N a 


1.3 (1.0,1.4) 


16.6 (15.9,17.2) 


1.0 (0.9,1.4) 


(1.4,2.5) 


W c /N s 


1.5 (1.2,1.5) 


15.8 (15.7,17.0) 


1.3 (1.3,1.4) 


(2.1,3.0) 


2Wc/N s 


1.1 (1.0,1.4) 


17.2 (16.2,17.4) 


1.7 (1.2,1.8) 


(1.7,2.2) 


published 


1.3 


16.54 


1.57 





TABLE I: Critical parameters from finite-size scaling. For 
com par sion, previously published results from Ref. [35j and 
Ref. |3ll are shown. In the first column, W c refers to the ac- 
cepted value of 16.5. Numbers in parenthesis are estimated 
bounds on parameters, and are based on the parameter re- 
gions Xrcd < 1- 



comparing 7 = 2W C /N S with 7 = W c /N s , we see that 
the results for larger 7 generally have larger uncertain- 
ties, but do not appear to be systematically shifted to- 
wards or away from their true values. It thus seems likely 
that one could obtain accurate values of W c and d 2 for 
larger values of 7 provided one can study systems that 
are large enough to keep the uncertainties to a reason- 
able size. As mentioned above, one has the additional 
requirement that 7 and A both have the same depen- 
dence on W in order obtain accurate values for v. This 
can be achieved, for example, by taking 7 oc l/po(W)N s , 
where po(W) is the ensemble- averaged density of states 
calculated for each strength of disorder. 



The next step is to obtain the critical exponent v, 
which is done by fitting a power law to 



dY 2 
dW 



-l/u 



no). 



(20) 



w=w c 



In Fig. HJc), we show dYijdW at W = W c , along with 
power law fits to the data. The three curves correspond 
to the best-fit, minimal, and maximal values of W c and d 2 
shown in Table |TJ The fitted exponents give three values 
of v for each 7, and are shown in the fourth column of 
Table |U We note that v is systematically underestimated 
for 7 tx W/N s , but is closer to the correct answer for 
7 oc W c /N s . One of the conclusions of this work is that 
v is more sensitive to the W^-dependence of 7 than either 
W c or d 2 . This follows directly from the derivative with 
respect to W in Eq. (|2"0)) , and means that obtaining an 
accurate value for v depends on establishing an accurate 
relationship between A and W. 

In summary, we have shown that the GIPR can distin- 
guish between localized and extended states, and more- 
over that it is possible to extract critical parameters from 
a scaling analysis of the GIPR. The main issue which 
arises is how the broadening 7 influences the results. In 
Table [Q comparing 7 = 2W/N S with 7 = W/N s and 



IV. CONCLUSIONS 

We have tested the usefulness of a generalized inverse 
participation ratio as a measure of Anderson localization 
by benchmarking it against the well-studied case of a dis- 
ordered three-dimensional tight binding lattice. Because 
the generalized inverse participation ratio depends on the 
local density of states, and not the single particle wave- 
functions, it is potentially useful for studying interacting 
systems where single particle wavefunctions are not de- 
fined. We have found that it is possible to extract critical 
parameters for the Anderson MIT, and have shown that 
finite size effects are not an impediment if the spectral 
broadening 7 used to calculate the local density of states 
is of the same order as the level spacing A. 
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